Measurement of the complete nuclide production and kinetic energies of the system 

136 Xe + hydrogen at 1 GeV per nucleon. 
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We present an extensive overview of production cross sections and kinetic energies for the complete 
set of nuclides formed in the spallation of 136 Xe by protons at the incident energy of 1 GeV per 
nucleon. The measurement was performed in inverse kinematics at the FRagment Separator (GSI, 
Darmstadt). Slightly below the Businaro-Gallone point, 136 Xe is the stable nuclide with the largest 
neutron excess. The kinematic data and cross sections collected in this work for the full nuclide 
production are a general benchmark for modelling the spallation process in a neutron-rich nuclear 
system, where fission is characterised by predominantly mass-asymmetric splits. 



I. INTRODUCTION 



In recent years, a vast experimental campaign has been 
dedicated to the measurement of spallation reactions at 
relativistic energies at the FRagment Separator [l| (GSI, 
Darmstadt). The installation of a target of liquid hy- 
drogen or deuterium 0] and the achromatic magnetic 
spectrometer [j| adapted to inverse-kinematics experi- 
ments were the tools for collecting high-resolution mo- 
mentum measurements and extracting production cross 
sections for each residue, identified in mass and nu- 
clear charge. Several systems, either favoured for un- 
veiling new physical aspects or directly relevant for ap- 
plications were studied. In spallation reactions, a large 
amount of the cross section which does not result in fis- 
sion fragments aliments the production of heavy nuclides; 
these heavy residues are, on average and almost inde- 
pendently of the neutron enrichment of the projectile, 
less neutron-rich than beta-stable nuclides. This prop- 
erty makes neutron-rich nuclei an ideal spallation tar- 
get for conceiving high- intensity neutron sources. Sys- 
tems like 197 Au (800 AMeV) +P II H, 2 ° 8 Pb(500AMeV) +P 
0,0, 208 Pb (1 ^ G eVl + P @,iLi 08 Pb (1 ^ G eV)+rf E3, 
JffUnA GeV) + P m E EE 13 El, 238 U(l A GeV) + d 

[la Il7|. were measured to study sequential evaporation, 
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and its intricated competition with fission. The mea- 
surement of 56 Fe+p at various energies focused 
the attention on systems below the Businaro-Gallone 
point [13, HH, in the fissility region where the saddle point 
becomes unstable towards asymmetric splits, and revived 
the discussion on the intermediate-mass-fragment forma- 
tion in spallation. 

Already twent y y ears ago, experiments in direct kine- 
matics [22l . [23l . I24I . I25L l2a | focused on the contribu- 
tion of this process in spallation; the corresponding phe- 
nomenological discussions [13, 0, 0, E3| suggested in- 
terpretations beyond the general fission-evaporation pic- 
ture (3ll [12] and in line with the onset of multifrag- 
mentation (reviews on this process can be found in 
Ref. A specific analysis of kinematic features con- 

nected to the intermediate-mass fragment formation in 
the 56 Fe( ly ^GeV) +P system [18| a nd in the much heav- 
ier system 238 U(iAGeV) + P uM pointed out the dif- 
ficulty of connecting this process exclusively to fission. 
In the former system, the presence of multifragmenta- 
tion was proposed as a relevant contribution. In the 
latter, the intermediate-mass-fragment formation was in- 
terpreted as fission events characterised by asymmetric 
splits, although surprisingly high fission velocities were 
observed. The modelling of the spallation reaction at 
relativistic incident energy depends largely on the degree 
of understanding of such a process. 

In this respect, the measurement of an intermedi- 
ate system was required to complete the survey on 
intermediate-mass-fragment formation. A new experi- 
ment was dedicated to the measurement of the com- 
plete residue production and the kinematic of the re- 
action 136 Xe(i a GeV) + P ■ 136 Xe is the stable nuclide 
with the largest neutron excess N — Z, with a fissil- 
ity below the Businaro-Gallone point. It is therefore 
best-suited for studying simultaneously the process of 
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intermediate-mass-fragment formation over a large range 
of light masses and the production of heavy evaporation 
residues from a system which has a neutron-to-proton ra- 
tio N/Z close to lead on the one hand, and in a fissility 
region where symmetric splits have low probability on 
the other hand. 

After a description of the experimental procedure, we 
present the measured cross sections for the production 
of fully identified nuclides formed in the reaction and 
the measured kinetic energies as a function of the mass 
of the residues. The compilation of cross sections and 
kinetic energies collected in the present work is the first 
to contain a large experimental survey on intermediate- 
mass-fragment production which extends to the heavy- 
residue production. 



TABLE I: List of layers interposed along the beam line in 
the target area and in the dispersive-plane region. 

Material Thickness [mg/cm 2 ] 



Target area: 
Vacuum window 
Beam-current monitor 
Mylar foils 

Front target windows 
Liquid hydrogen 
Rear target windows 
Mylar foils 
Dispersive plane: 
Scintillator 
Degrader (wedges) 



Ti 4.5 
Ti 13.5 
C5H4O2 4.15 
Ti 18.15 
H 2 87.3 
Ti 18.15 
C5H4O2 4.15 a 

C9H10 475.45 (C) + 44.025 (H) 
Al 816. 6 b 



"Mylar is coated with 0.1 mg/cm 2 of aluminium 
'For the reference trajectory. 



II. 



EXPERIMENTAL PROCEDURE 



The fragments were produced in inverse kinematics by 
directing a primary beam of 136 Xe at 1 A GeV on a target 
of liquid hydrogen contained in a cryostat with thin tita- 
nium windows. The projectile residues were then anal- 
ysed in-flight, using the inclusive measurement of their 
momenta along the beam axis. The placement of the de- 
tectors in the FRagment Separator for this experiment 
is sketched in Fig. [TJ In the target area a beam-current 
monitor was placed to measure the primary-beam inten- 
sity. The positions where the trajectories of the frag- 
ments intersect the dispersive focal plane (sdfp) and 
the terminal focal plane (xtfp) were registered by two 
scintillation detectors. Their combined signals provided 
in addition the measurement of the time of flight. The 
fragments of 136 Xe, which at relativistic incident energies 
are fully stripped with high probability, were identified 
in nuclear charge Z by two ionisation chambers placed in 
front of the terminal focal plane. 

The primary beam did not only interact with the hy- 
drogen contained in the cryostat, but also with the cryo- 
stat itself, wrapped into insulating Mylar foils, and with 
the other layers of matter present in the target area, like 
the beam-current monitor and the accelerator vacuum 
window. To measure the contribution of non-hydrogen 
target nuclei, the experiment was repeated under equal 
conditions, by replacing the liquid hydrogen target with 
an identical empty sample. Table U lists the composition 
and the thicknesses of all layers of matter placed in the 
beam line during the experiment. 



High-resolution achromatic mode 



with opposite dispersion factors. The spectrometer has 
maximum dispersion in the centre of the dispersive fo- 
cal plane (xdfp) and was set to be achromatic. More 
precisely, xdfp is related to the momentum deviations 
5p = (p — po)/po with respect to the reference trajectory 
Po, in the first dispersive section, and 5p', in the second 
dispersive section, by 



D5p + gxi F p = .tdfp = D'Sp' + g'xTFP, 



(1) 



where Xifp and xtfp are the displacements in the initial 
and terminal focal plane, respectively, from a reference 
trajectory which, by convenience, was chosen to intersect 
all focal planes in their centres. D and D' are dispersion 
constants; g and g' are the magnification factors mea- 
sured when moving from the extremes of the beam-line 
towards the dispersive focal plane. The three optical pa- 
rameters D, D' and g' were measured with uncertainties 
of 0.7%, 2.7% and 1%, respectively, in an initial cali- 
bration run and later kept fixed to constant values for 
the whole experiment. We did not need to measure g 
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In essence, the magnetic configuration of the spectrom- 
eter is based on four bending sections delimited by an 
initial focal plane, where the target is placed, and a ter- 
minal focal plane. The first pair of bending magnets and 
the second pair, with average magnetic fields B and B' , 
respectively, constitute two portions of the spectrometer 



FIG. 1: (Color online) Layout of the FRagment Separator. 
The positions of the four dipole magnets, the focal planes and 
the main detectors are shown in a horizontal plane view, in 
scale along the beam direction. The limited angular accep- 
tance of ~ 15 mrad in the laboratory frame is indicated. 
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as the beam hits the target in the centre of the initial 
focal plane. As eq. (TJ shows, the momentum deviation 
of the fragments is completely denned by the displace- 
ments xdfp and xtfp, which were measured by placing 
a scintillator detector in each of the corresponding focal 
planes. By substituting the momentum deviation Sp by 
the magnetic-rigidity deviation 8 {Bp) = {Bp—Bp )/Bp 
with respect to the reference trajectory Bp , eq. leads 
to the two equations that govern the data analysis: 

BP = B.p (l + ^), (2) 

{BP)' = B'.p' (l + X ™-J X ™ ). (3) 

The curvature radii po and p' were kept fixed for the 
whole experimental run. In order to scan the full distri- 
bution of magnetic rigidities of the fragments, the mag- 
netic fields B and B' had to be changed several times 
due to the limited acceptance in magnetic rigidity, which 
selects ranges of around 3% for each individual magnetic 
setting {B,B') in the dispersive focal plane. In order to 
keep all the optical parameters strictly unchanged, the 
magnetic fields of the ensemble of magnets were scaled 
by equal factors for the first and the second dispersive sec- 
tion, respectively. These two factors should differ slightly 
in order to keep the range of elements selected by the 
spectrometer fixed for a set of measurements. 



B. Separation of fragments 

As the time of flight could be measured between the 
dispersive and the terminal focal plane, we could asso- 
ciate a mass-to-charge ratio A/Z to each fragment having 
magnetic rigidity {Bp)' in the second dispersive section: 

± = l c M (4 ) 

Z c mo + 8m (3j 

where c is the velocity of light, e the electron charge mag- 
nitude, m the nuclear mass unit, 8m — &M/A the mass 
excess per nucleon, and (3-y the relativistic factor, which 
is determined in the laboratory frame in longitudinal di- 
rection. 

Due to the limited acceptance in magnetic rigidity, the 
scanning of the whole A/Z distribution of fragments re- 
quired several magnetic-field scalings. When thick layers 
of matter are inserted between the two dispersive sec- 
tions, the fragments lose part of their kinetic energy as 
a quadratic function of the charge, and their magnetic 
rigidities change. This property can be exploited[3J|: it 
imposes a charge-selection in the second dispersive sec- 
tion which can be employed to measure restricted groups 
of elements. For this purpose, we used an aluminium 
layer (degrader) of 816.6 mg/cm 2 in the beam line and 
selected three bands of nuclides centred around silver, 
zinc and aluminium, according to the expectation that 



the production yields do not vary largely within each sin- 
gle band. The aluminium degrader was shaped in order 
not to perturb the achromatic mode. 



III. ANALYSIS 

The first step of the analysis consisted of identifying 
each fragment in nuclear charge Z by calibrating the 
ionisation chambers and in mass A from the magnetic 
rigidity and the time-of-flight measurement of their mo- 
mentum in the second dispersive section of the spectrom- 
eter. Afterwards, from the knowledge of the A/Z ratio 
and the magnetic rigidity measured in the first dispersive 
section of the spectrometer, the momentum was deduced 
a second time with higher precision than from the time- 
of-flight measurement. Hence, a high-resolution distri- 
bution of longitudinal recoil velocities in the projectile 
frame uj? was associated to each reaction product iden- 
tified in nuclear mass and charge. These distributions 
were normalised to the beam dose per target thickness, 
and the parasitic effect of the reactions in layers differing 
from hydrogen was removed, so that the integral of the 
distribution was equal to the measured production yield 
/ for each nuclide. The analysis then focused on the 
shape of these spectra which, without a dedicated analy- 
sis procedure, were still not directly suited for extracting 
the physical quantities related to the kinematics and the 
production of the reaction. In order to reconstruct the 
full distribution of emission velocities, independent of the 
experimental conditions, the geometry of the angular ac- 
ceptance of the spectrometer was accurately accounted 
for. This reconstruction was necessary for extracting the 
production yield of each nuclide as well as the kinetic- 
energy distributions. Finally, the formation cross section 
a for each primary reaction product was extracted from 
the production yield by accounting for the secondary re- 
actions induced by the primary fragments in the hydro- 
gen target and in the dispersive plane, where a scintillator 
detector and a degrader were placed. 



A. Nuclide identification 

The relativistic factor f3 in eq. {!]) is the ratio (i = l/ct, 
where I is the path length, which is given by £q = 36 m for 
a fragment centred at all focal planes, and t is the actual 
time of flight, t could not be deduced directly from the 
measured time of flight t' due to non-linear effects of the 
light-propagation time in the scintillating detectors (re- 
sulting into quadratic terms in the path xp,pp and xtfp) 
and an amplitude dependence on the charge of the signal. 
In order to extract the actual values of I and t, a set of 
eight coefficients fcj, was introduced in the form: 

£ = £ + kiXuFP + fe^TFP, (5) 
t = k 3 + k^' + k 5 xl, F p + k 6 x^ FP + k 7 e~ k8z2 . (6) 
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FIG. 2: (Color online) Ensemble of all events, identified in 
nuclear charge and mass. The overlapping bands a, b and c 
correspond to the three groups of magnetic settings for nu- 
clide distributions centred around 120 Ag, 69 Zn and 24 Al, re- 
spectively. The band c, corresponding to light nuclides, is 
enlarged to emphasise the quality of the isotopic resolution. 



The terms fc, were deduced by numerical optimisation 
and used for the whole data analysis. 

The full nuclide identification was obtained from 
eq. ([3]) and eq. ([4]). The raw data are shown in Fig. [21 
where all the events collected in the experiment are shown 
as a nuclide identification plot. 



B. Beam dose per target thickness 



where ab is a coefficient to convert the secondary-electron 
current produced by the primary beam into the number 
of projectiles impinging on the target, fb(i) is a measure- 
ment of the primary-beam current by the beam-current 
monitor [35[, t(i) is the measured relative dead time of 
the data acquisition system and cvh 2 is the number of 
nuclei per area of the liquid hydrogen. 

The coefficient £tb was deduced in a dedicated measure- 
ment: the beam-current monitor was calibrated in com- 
parison with a scintillation detector following the method 
described in ref. [3f3]. The spill structure digitised by the 
two detectors is shown in Fig. 03 After accurately sub- 
tracting the offset of the beam-current monitor and inte- 
grating the recorded counts over each spill, a quadratic 
dependence of the particle counting as a function of the 
beam intensity was obtained due to the saturation of the 
scintillation detector, as shown in the inset of Fig. 03 
The coefficient ah was deduced as the initial slope of the 
function, with an uncertainty of 1%. In a previous exper- 
iment [37| a slight dependence of this coefficient with the 
position of the beam spot on the target (which has a fluc- 
tuation of the order of one millimetre) was estimated to 
introduce an additional uncertainty of around 5%, which 
we apply to the present data. 
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FIG. 3: (Color online) Beam-current monitor calibration. 
Superposition of counts in the beam-monitor (coloured spec- 
trum, axis label on the left) and counts in the scintillator 
(white-filled spectrum, axis label on the right). The beam- 
monitor spectrum is multiplied by the parameter Ob, that 
coincides with the calibration slope (solid line) shown in the 
inset. 



The number of events registered in different settings 
of the separator magnets were normalised to the same 
beam dose in order to have consistent weights. The nor- 
malised counts N(i), registered for an individual experi- 
mental run i, determined by a specific magnetic setting, 
are obtained dividing the number of events n(i) by four 
coefficients: 



ab/b(i)[l - r(i)]aa 2 ' 



C. Longitudinal recoil velocities 

Once a fragment is identified in mass and charge, its ve- 
locity is directly obtained from the magnetic rigidity. As 
no further layer of matter was present behind the target 
in the whole first dispersive section of the spectrometer, 
the relativistic factor (3j in the laboratory frame in lon- 
gitudinal direction could be deduced more precisely from 
the magnetic rigidity Bp, defined in eq. ((2]), rather than 
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from {Bp)' , defined in eq. ([3]), so that 

1 e Z 



= Bp 



c mo + 5m A 



(8) 



To change from the (/?7)| factor in the laboratory 
frame to the longitudinal velocity ik in the beam frame, 
the energy loss in the target was taken into account. In 
particular, we assumed that the average position where 
the collision occurs corresponds to half of the total thick- 
ness of the ensemble of layers present in the target area 
(listed in table HJ ■ We defined the beam frame in cor- 
respondence with the velocity of the projectile at this 
position. The high precision in deducing recoil velocities 
for individual reaction products is given by the measure- 
ment of Bp (or (/37)]j J ), which has a relative uncertainty 
of only 5-1CT 4 (FWHM). However, to access the reaction 
kinematics, the details of the shape of the longitudinal- 
velocity distribution d-/V(i>|| )/duf| had to be examined for 
each nuclide. 

In section [I] we explained that the spectrometer has a 
limited acceptance in momentum; for this reason the full 
longitudinal momentum distribution for a single nuclide 
was constructed by composing several measurements ef- 
fectuated with different magnetic settings. The spec- 
trometer has also a limited acceptance for the emission 
angles so that it can be demonstrated that the measured 
spectra for intermediate-mass fragments are close to in- 
variant cross sections This problematic is revisited 
in detail in section lTlIDI Another technical difficulty was 
the contribution to the spectrum diV(vjj ) )/di;jj ) of reac- 
tions occurring in any layer of the target area other than 
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FIG. 4: (Color online) Subtraction of the contribution of 
non-hydrogen target nuclei (label "d", for "dummy") to obtain 
the measured velocity distributions for n B, 48 Ti and 114 Sn. 
The left spectra represent the full measured contribution to 
the velocity distribution associated to the hydrogen target 
(colour-filled histograms) and to non-hydrogen target nuclei 
(black-filled histograms), while the right spectra represent the 
deduced hydrogen contribution alone. The contribution from 
the windows is large for light fragments ( n B, 48 Ti) and neg- 
ligible for heavy residues ( 114 Sn). 



liquid hydrogen. To solve this problem, all experimen- 
tal steps were repeated under equal conditions with an 
empty-target sample in order to deduce the longitudinal- 
velocity distribution diV'(ujj ) )/d?;jj ) , related to parasitic re- 
actions in the target area in all layers other than hydro- 
gen. The measured distribution of yields d/(vj| )/di>j| was 
directly obtained as the difference between the two distri- 
butions dN(v^ )/duy and dN'(v^)/dv^, both normalised 
to the number of nuclei per area of liquid hydrogen cvh 2 • 
It should be remarked that, while the parasitic layers 
of matter have been designed to be as thin as possible so 
as to maximise the relative production in hydrogen, they 
still induce, on average, more violent reactions, resulting 
in the production of intermediate-mass fragments with 
large yields. As a consequence, in the intermediate mass 
range, the parasite contribution can exceed 50%. This 
fact, illustrated in Fig. |4j required the same accuracy 
for the measurement of the parasite reaction as for the 
production in the full target. In the case of nuclides 
produced by reactions of charge exchange, the parasite 
contribution was smaller than 4% for Cs isotopes and 
smaller than 5% for Ba isotopes. 



D. Reconstruction of the angular acceptance 

By employing the method exposed so far, for each 
individual nuclide we measured a longitudinal-velocity 
spectrum dI(vYA/dv^. This is the distribution of the 

longitudinal-velocity component vj? in the projectile 
frame corresponding to the portion of the velocity-space 
distribution da/dv h selected by the angular acceptance. 
The velocity-space distribution could be reconstructed by 
a deconvolution procedure [ll|, directly from the mea- 
sured spectrum dl(v^)/dv^. Fig. [5] presents in the left 
column the measured longitudinal-velocity spectra for 
three intermediate-mass fragments. The corresponding 
reconstructed velocity-space distributions are shown in 
the right column by planar cuts along the beam axis. 
The boundary lines on the planar cuts indicate the angu- 
lar acceptance of the spectrometer. The probability that 
a fragment is emitted within these boundaries is defined 
as the transmission probability, which is equal to unity 
for the heaviest residues and drops to smaller values for 
light fragments. 

For light elements up to silicon, the reconstruction of 
the velocity-space distribution was effectuated in rela- 
tion with two components which could be disentangled 
in the measured longitudinal-velocity spectra, as illus- 
trated in Fig. [H In general, we observe that the shape 
of the longitudinal- velocity spectra for the intermediate- 
mass fragments evolves between two extreme patterns, 
characterised by a concave and by a convex centre, re- 
spectively. Physically, these two modes reflect the action 
of the Coulomb field, which acts in different ways depend- 
ing on the decay process. The velocity-space distribution 
of the concave mode is an isotropic shell in the frame of 
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the emitting source; this is the effect of the Coulomb re- 
pulsion in decay processes with low fragment multiplic- 
ity and high fragment-size asymmetry like, for instance, 
asymmetric splits in fission or multifragmentation events 
characterised by a low multiplicity of fragments. The 
convex mode can be associated with two very different 
processes: either residues of sequential evaporation, or a 
much more violent process of multifragmentation in sev- 
eral pieces of comparable size. The asymmetry of the 
concave component can be fully described as a trivial 
effect of the angular acceptance (the integral of the for- 
ward peak is 14% larger that the integral of the backward 
peak for n B and 12% for 31 Si). On the other hand, we 
attribute the asymmetry of the convex component mainly 
to the mixing of emission processes associated with dif- 
ferent sources: we assume in fact that a large range of 
excitation energies involved in the formation of a given 
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FIG. 5: (Color online) Left column. Measured velocity dis- 
tributions for three nuclides: n B, 20 F, 31 Si. All spectra are 
normalised to the same integral value. The spectra are fitted 
(green curves) in order to determine the concave (red curves) 
and the convex (blue curves) components; see text. Right col- 
umn. Planar cuts along the beam axis of the reconstructed 
velocity-space distribution dcr/d?; b in the beam frame. All 
distributions, reconstructed from the corresponding measured 
velocity distributions, are normalised to the same integral 
value and described by a logarithmic evolution of the colour, 
where the span from blue to pale red corresponds to a factor 
of two in the intensity. The two lines indicate the boundaries 
of the angular acceptance, inside of which the fragments could 
be measured. 



intermediate-mass fragment would be reflected in a large 
range in momentum transfer, which is a quantity related 
to the violence of the reaction j38|. 

The deconvolution procedure is performed assuming 
that the velocity-space distribution is composed of sev- 
eral emission processes, which are symmetric with respect 
to the longitudinal axis; we associate the concave mode 
to one isotropic source and we describe the asymmetric 
shape of the convex component by several sources with 
the same Gaussian shape and different weights. For de- 
tails see ref . [3§| . The distribution of the emitting sources 
has to be deduced by an optimisation procedure. The 
hypothesis of isotropic emission from a given source is 
adapted to relativistic collisions |l8| and would not be 
valid any more at Fermi energies. As shown in Fig.[6j the 
transmission probability obtained by the deconvolution 
procedure is very different for the two kinematical modes. 
Also the uncertainties differ remarkably: the large er- 
ror bars for the convex mode reflect the uncertainty in 
deducing the distribution of emitting sources and, more 
generally, the difficulty in determining the width of the 
distribution of the transverse velocity associated to this 
mode. In contrast to this, the concave mode is well- 
controlled even at much lower transmission. In Fig.[6]the 
mass evolution of the transmission probability for inter- 
mediate mass fragments is shown: it is the composition 
of the transmission probability for the concave and con- 
vex modes and is labelled as "multiple-source" approach, 
to indicate that the deconvolution was performed under 
the assumption that several sources were involved. 

For elements above silicon the two kinematic modes 
can not be disentangled any more because the measured 
integral of the concave mode decreases with respect to 
the convex mode and the spacing between the two humps 
of the concave mode reduces. Besides, the method de- 
scribed so far is necessary if the measured spectra are 
characterised by a large dispersion of emitting sources, 
and if, in particular, the convex mode shows a sizable 
asymmetry and is globally displaced with respect to the 
concave mode. In this respect, above silicon, this method 
is much less justified because the whole shape of the 
measured spectra tends to be symmetric and Gaussian 
(this trend is visible in Fig. [5]). We assumed in this case 
that all sources alimenting the emission of a given frag- 
ment coincide and that the kinematics is fully isotropic. 
This simplification was applied in the same manner in the 
analysis of light fragments produced in the spallation of 
iron where the source velocities of the two kinemati- 
cal modes contributing to the production of fragments of 
a given mass were about the same. The mass evolution of 
the transmission probability calculated with this simple 
assumption is shown in Fig. [6] and is labelled as "unique- 
source" approach, to indicate that the deconvolution was 
performed with the assumption that the whole kinemat- 
ics could reduce to one single source. The smooth func- 
tion is the result of a fit and the uncertainty contains the 
scattering around the function. Such uncertainty reflects 
the quality of the experimental velocity spectra, which 
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Transmission probability 




A 



FIG. 6: (Color online) Evolution of the transmission prob- 
ability as a function of the mass of the fragments, deduced 
assuming their emission geometry to be either related to one 
unique source or related to more sources. In the latter case, 
the kinematics is assumed to be the composition of the con- 
cave and the convex mode. The transmission evaluated sepa- 
rately for the two kinematical modes is also shown. 



are affected by low statistics in the region of the lowest 
production yields (around half the mass of the projectile); 
it reflects also the numerical difficulties in the deconvolu- 
tion when the transmission probability approaches unity 
(for A w 80): in that case, the shape of the angular ac- 
ceptance affects the results strongly. Above A « 100 all 
fragments match the acceptance and the calculation is 
trivial. For comparison, the transmission probability for 
fully isotropic emission is also calculated for the light- 
est fragments and compared with the "multiple-source" 
approach in Fig. [6l The difference between the two ap- 
proaches is up to 30%. 



E. Extraction of primary production cross sections 



To extract the production cross section from the pro- 
duction yields, we had to correct for the occurrence of 
secondary reactions in the layers of matter present along 
the beam-line. Secondary-reaction products formed in 
the dispersive plane have different consequences on the 
measurement than those formed in the target area. 

A secondary-reaction product formed in the dispersive 
plane deviates from the trajectory related to the mag- 
netic rigidity of the corresponding mother nucleus, it is 
spatially separated by the ion optics and, with high prob- 
ability, not transmitted. We corrected for the loss of the 
primary production by calculating the attenuation of the 
beam of fragments when traversing the scintillator de- 
tector and the degrader. The probability for a nucleus 
(j4oj Zq) to traverse a layer of matter of x atoms per area 
without interacting is equal to Vq = exp[— a Q ■ x], and 
depends on the total reaction cross sections ctq. The to- 



tal reaction cross sections were calculated according to 
the model of Karol [5(J modified by Brohm [4l|, with an 
uncertainty of 5%. The correction to apply as a multi- 
plicative factor to the measured yields is equal to 1/Vo 
and varies from around 1.05 for the lightest fragments 
to 1.15 for the heaviest residues; it is illustrated as a 
function of the mass number of the measured residue in 
Fig. Eta). 

When secondary-reaction products are formed in the 
target area (liquid-hydrogen target), suppressing the sec- 
ondary fragments by the magnetic spectrometer is not 
possible because they are produced before entering both 
dispersive sections: at the same time and with no dis- 
tinction we measure a slightly reduced distribution of 
primary reaction products together with a distribution of 
secondary reaction products. No direct experimental ob- 
servables can be related to the loss and gain of production 
yields due to secondary reactions in the target area. Even 
though formally we can establish exact relations between 
the primary- and secondary- reaction fragments, these re- 
lations require the knowledge of the nuclide-production 
cross sections in the reaction between primary fragments 
and the target. The correction is therefore dependent on 
the reaction model we apply. 

The nuclide-production cross sections in the sec- 
ondary reactions were calculated by coupling the mod- 
els BURST and ABLA [II El, 0] ; the former is a 
parametrisation of the intra-nuclear cascade, the latter 
models the decay of hot fragments. To calculate the cor- 
rection factor we used the same method of Ref . [HI , em- 
ployed in the analysis of a recent experiment on the spal- 
lation of 238 U [l3|. The method evaluates the secondary- 
production yields in a thick target in inverse kinematics 
at relativistic energies. In the analysis of the spallation of 
238 U the method considered decays by evaporation and 
fission; in the present analysis, we employed this method 
in a reduced form, excluding fission, even though we ex- 
pect its contributions in the production of intermediate- 
mass fragments in the spallation of xenon. The modelling 
of these particular channels are not well controlled at this 
incident energy due to a lack of experimental data. For 
this reason, the correction factor was set equal to unity 
for light fragments. Since heavy fragments are mostly 
produced by evaporation, their correction factor is not 
much affected by this limitation. The result of the cal- 
culation is shown in Fig. Ulb) on a nuclide chart. The 
correction factor (multiplicative coefficient applied to the 
measured production yields) accounting for secondary re- 
action in the target area is set to 1 for the lightest frag- 
ments up to aluminium. It then decreases gradually from 
unity down to about 0.8 for the isotopes of nickel; it in- 
creases monotonically for heavier elements and in the re- 
gion of zirconium it exceeds unity, up to the largest value 
1.06 attained in proximity of 136 Xe . 

The global effect of the secondary reactions in the tar- 
get is to reduce the yields of the heaviest elements in 
favour of those fragments having around half the mass of 
the projectile, which are also the least produced. More in 
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detail, due to the tendency of evaporation residues to mi- 
grate towards the neutron-deficient side of the chart of 
the nuclides, secondary reactions substantially increase 
the yields of the less neutron- rich fragments. This ten- 
dency is illustrated in Fig.[7]Jb) and in the corresponding 
inset, where the evolution of the correction factor is stud- 
ied for four elements as a function of the neutron num- 
ber. In the case of Cs and Ba the loss of cross section 
due to secondary spallation reactions producing lighter 
residues is compensated by the gain of cross section due 
to secondary reactions of charge exchange. The uncer- 
tainty on the correction factor depends on two quantities. 
First, it depends on the uncertainty of the calculation of 
the total reaction cross sections which, also in this case, 
were obtained from the model of Karol [13] modified by 
Brohm [4l|, with an uncertainty of 5%. Second, it de- 
pends on the uncertainty that we attribute to the model 
calculation of the nuclide-production cross sections in the 
secondary reactions: according to the simulation of mea- 



Correction for secondary reactions in the dispersive plane 
1.2 r 




2 8 J Rejected 



FIG. 7: (Color online) (a). Correction factor for secondary 
reactions in the dispersive plane (attenuation in the scintil- 
lator and the degrader) as a function of the mass number of 
the residues. The uncertainty is indicated by the coloured 
band. (b). Correction factor for secondary reactions in the 
target area (liquid hydrogen) on a nuclide chart; in the inset 
it is illustrated for four elements as a function of the neutron 
number. 



sured spallation data (see e.g. Ref. [14]), we estimated 
that this uncertainty is as large as 20%. This second 
contribution most strongly affects the nuclides which are 
largely produced by secondary reactions and the uncer- 
tainty is larger for nuclides associated to smaller correc- 
tion factors, as illustrated in the inset of Fig. [3(b) . 

The experimental method turns out to be limited for 
the measurement of neutron-deficient residues as their 
yields are completely dominated by secondary reactions. 
Yields alimented for more than 50% by secondary reac- 
tions were rejected. 



IV. EXPERIMENTAL RESULTS 

The average kinetic energies were deduced from the 
velocity-space distributions reconstructed da/dv h . From 
the same unfolding procedure, as discussed in sec- 
tion [ThEI the nuclide cross sections were obtained from 
integrating the reconstructed velocity spectra and cor- 
recting for secondary reactions, so that the production 
cross sections were deduced as: 

° = I ■ fir ■ /tar ' /dP, (9) 

where / tr is the transmission probability, f tar the cor- 
rection for secondary reactions in the target and /dp the 
correction for secondary reactions in the dispersive plane. 



A. Kinetic energies 

A general survey of the average kinetic energy im- 
parted to the spallation residues is presented as a func- 
tion of the mass number in Fig. [8l and the data are 
listed in table HH It results from the decay kinematics 
as well as from the displacement of the emitting source 
in the beam frame as a consequence of the collision. The 
latter contribution is shown in a separate spectrum in 
Fig. [H The average kinetic energies were deduced from 
the velocity-space distribution reconstructed according 
to the "multiple-source" prescription described in sec- 
tion HUE] up to A = 30. The error bars reflect the un- 
certainty in the fit of the measured longitudinal- velocity 
distribution, the uncertainty in the unfolding procedure 
and the statistics of the measurement. The enlarging of 
the error bars in approaching A = 30 reflects the increas- 
ing difficulty in applying the "multiple-source" prescrip- 
tion to the measured longitudinal-velocity distribution 
which, becoming gradually closer to a Gaussian distribu- 
tion, progressively reduces the information for estimating 
the asymmetry of the emission kinematics. For masses 
larger than A = 30 we switched to the "unique-source" 
prescription as the measured velocity spectra are in gen- 
eral Gaussian-like. The mismatch at A = 30 indicates 
that the kinematics is still not isotropic in this mass re- 
gion and the two approaches are not yet equivalent. 



TABLE II: Kinetic energies in the projectile frame averaged 
over isobaric chains. The uncertainty, indicated in paren- 
thesis, includes both statistical and systematical errors and 
applies to the last decimal digits. 
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For the heaviest fragments, which are characterised 
by measured longitudinal velocity spectra with Gaus- 
sian shape, of width a v and mean value (uj?), the kinetic 

energy is given by E k = (l/2)m A((wj^) 2 + 3cr 2 ). The 
heaviest fragments are qualitatively consistent with the 
systematics of Morrissey which was deduced from 
residues of sequential evaporation in very peripheral ion- 
ion collisions. This systematics correlates the mean value 
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FIG. 8: (Color online) The distribution of mean kinetic en- 
ergies in the projectile frame averaged over isobaric chains is 
represented by the upper diagram. The lower spectrum indi- 
cates the contribution to the mean kinetic energy that is at- 
tributed to the mean momentum transfer during the impact. 
The difference between the two spectra is the kinetic energy 
gained in the decay. The error bars include both statistical 
and systematical errors. A vertical line separates the por- 
tion of the upper spectrum related to fragments lighter than 
A — 30, evaluated according to the "multiple-source" prescrip- 
tion, from the rest of the spectrum, evaluated according to the 
"unique-source" prescription. The diagram is divided in two 
parts corresponding to half the mass of the projectile. In the 
heavier-mass portion both experimental spectra are compared 
with the Morrissey systematics [38[] (dashed line). 

and the width of the momentum-transfer distribution in 
the projectile frame to the mass loss and to the square 
root of the mass loss, respectively. 

B. Production cross sections 

The experimental production cross sections are col- 
lected in table Hill and IIVI and shown as isotopic distribu- 
tions in Fig. [9] and Fig. [TUl Despite an effort in measur- 
ing the cross sections down to the microbarn level, the 
measurement of very low cross sections was, for neutron- 
deficient nuclides, hindered by large feeding from sec- 
ondary reactions in the target. The uncertainty, rarely 
exceeding 15% for the lightest nuclides and below 10% 
on average for the heavy nuclides accounts for all statis- 
tical and systematical contributions: statistics, analysis 
and fitting of the velocity spectra, the procedure of ve- 
locity reconstruction (reflected in the uncertainty of the 
transmission probability) , calibration of the beam mon- 
itor, secondary reactions in the target and in the beam 
line. The formation cross sections were deduced for each 
measured nuclide, for the intermediate-mass fragments, 
for the heavy residues, and even for the production of 
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TABLE III: Isotopic cross sections for the production of elements ranging from Li to Tc. The uncertainty includes statistical 
and systematical errors. Where indicated in parenthesis, it applies to the last decimal digits. 
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0.1747(99) 


95 


n /^^t/o^\ 

0.647(34) 


■t A 1 /ir/i n \ 

14 1.65(12) 


32 


n nnni /i t\ 

0.0201(17) 


52 


n n i n 1 / 1 /"> \ 

0.0104(16) 


62 


n i o n / 1 \ 

0.182(16) 


74 


n nTn / 1 o \ 

0.270(18) 


87 


n nTno / A -i \ 

0.0708(41) 


96 


n 07i/io\ 

0.351(18) 


i 7 n 1 T/""/nT\ 

15 0.176(27) 


33 


n ii77/T7.\ 

0.1155(75) 




V 


63 


n i no / 1 n\ 

0.103(10) 


75 


pi i nn / 1 o \ 

0.199(13) 


88 


n no 7n/ 1 c\ 

0.0272(16) 


97 


n i to i / nn\ 

0.1731(90) 


N 


34 


n o o 7. / 1 o \ 

0.285(18) 


46 


0.00030(12) 


64 


0.0561(66) 


76 


n i nor / /"' t \ 

0.1025(67) 


89 


0.01024(68) 


98 


pi nTno/ a o\ 

0.0792(42) 


13 0.0241(28) 


35 


n n ^ n / i /"* \ 

0.240(16) 


47 


0.0099(14) 


65 


n aata/h \ 

0.0250(31) 


77 


n ataa/o d\ 

0.0500(34) 


90 


n nno/^i /ni \ 

0.00361(31) 


99 


n nooo/inN 

0.0338(19) 


1^ n H /I/ 1 / 00\ 

14 0.446(33) 


36 


n 1 rn/i n\ 

0.159(12) 


48 


n n 1 n n / o \ 

0.0492(63) 


66 


n nnn a ( i i \ 

0.0094(11) 


78 


n mno/i ^\ 

0.0198(14) 




Sr 


100 


n ni ooo/on\ 

0.01338(82) 


1 f 1 or/i o \ 

15 1.85(13) 


37 


n at ai f v p \ 

0.0591(56) 


49 


n i o t / 1 t\ 

0.137(17) 


67 


n n n nrn/ 

0.00352(48) 


79 


n nnTn/^ / p r\ \ 

0.00706(60) 


81 


n no /nn\ 

0.084(22) 




Mo 


i/ -1 n /i nn / o o '\ 

16 0.499(38) 


38 


r\ nn a a / o n \ 

0.0244(32) 


50 


n i Tn / aa\ 

0.172(20) 




Cu 


80 


n nnnn i / /i ^ \ 

0.00291(44) 


82 


0.360(46) 


90 


n p\p r\ f a 

0.269(47) 


it n o^c / cn\ 

17 0.276(50) 




CI 


51 


Pi 1 /II / i o\ 

0.141(18) 


60 


0.00845(71) 


81 


Pi P\P\~\ P\ A 1 rro\ 

0.00104(52) 


83 


n Tn n\ 

0.795(59) 


91 


i one/ pict \ 

1.205(96) 





34 


0.00737(61) 


52 


n nTi /i i \ 

0.071(11) 


61 


n Aror/ ho\ 

0.0585(43) 




Se 


84 


i n v o / /""tN 

1.058(67) 


92 


2.54(16) 


i v n no n / 1 o '\ 

15 0.039(13) 


35 


n nor"^/r*T\ 

0.0854(57) 


53 


Pi P\Af\0/*7G\ 

0.0403(78) 


62 


n ii~t/ii\ 

0.157(11) 


72 


n n/^r"o/nT\ 

0.0653(97) 


85 


i n ^ n / n\ 

1.049(60) 


93 


o i / i t\ 

3.16(17) 


1 p n p no / a n\ 

16 0.603(40) 


36 


n 1 o d / 1 n\ 

0.184(12) 


54 


0.0143(41) 


63 


n nn i / 1 \ 

0.221(16) 


73 


n nnn / 1 o \ 

0.229(18) 


86 


n onn//i/^\ 

0.822(46) 


94 


n n / i a \ 

2.62(14) 


IT Pi CO A 1 00\ 

17 0.534(38) 


37 


0.245(16) 


55 


Pi PiPi C il /oo\ 

0.0054(23) 


64 


0.194(14) 


74 


Pi A oo /qa\ 

0.422(30) 


87 


n toe/ on\ 

0.526(29) 


95 


2.05(11) 


io n p r p / 7 n\ 

18 0.656(50) 


38 


r\ inn/inN 

0.139(10) 




Cr 


65 


n inno/no\ 

0.1203(98) 


75 


n /inT/oi\ 

0.407(31) 


88 


n nTT / 1 r\ 

0.277(15) 


96 


i nnn//" 1 t\ 

1.299(67) 


in n n 7 /"> / o o \ 

19 0.256(33) 


39 


n nTT a / pd\ 

0.0774(68) 


49 


n nn,in/in\ 

0.0049(12) 


60 


0.0689(66) 


76 


n onn/no\ 

0.322(28) 


89 


n inoo/Ti\ 

0.1283(71) 


97 


n n i n / A t\ 

0.912(47) 


F 


40 


r\ n o f o / ^n\ 

0.0358(40) 


50 


n a d 1 f / n n \ 

0.0415(92) 


67 


n AOA1 / yl1 \ 

0.0321(41) 


77 


n nno / 1 n\ 

0.208(19) 


90 


n nr*in/nn\ 

0.0512(29) 


98 


n r'^7^/on\ 

0.574(30) 


it n nnnrir / Tn\ 

17 0.00925(79) 




Ar 


51 


n i o i / nT\ 

0.131(27) 


68 


n pfi r\p / -i p\ 

0.0126(16) 


78 


n 1 ta/i a\ 

0.150(10) 


91 


n pr\r\r\ ( i n\ 

0.0202(12) 


99 


n o n o / 1 /""■ \ 

0.303(16) 


18 0.095(10) 


36 


n nn 7 n / 1 1 \ 

0.0059(11) 


52 


pi i ne f ap\\ 

0.196(40) 


69 


0.00572(62) 


79 


pi noon /ci \ 

0.0820(51) 


92 


n nnTno/cr.cr\ 

0.00792(55) 


100 


0.1451(75) 


in n o o o / n o \ 

19 0.388(28) 


37 


n Plf P A f A A \ 

0.0564(44) 


53 


n i 1 n / o i \ 

0.149(31) 


70 


n nni or 1 / ,i n \ 

0.00136(42) 


80 


0.0362(23) 


93 


n nnnTn/o^N 

0.00270(34) 


101 


n n/^To/o/^N 

0.0678(36) 


nn n ^ n ^ / A A \ 

20 0.494(44) 


38 


n i o f / i n \ 

0.185(12) 


54 


0.096(20) 




Zn 


81 


n m ono/nn\ 

0.01303(92) 




Y 


102 


n non/" 1 /i t\ 

0.0306(17) 


n -i n /i n ^ /Tn\ 

21 0.404(70) 


39 


n ni o/i 7 N 

0.218(15) 


55 


0.0410(84) 


62 


n nn^n/io\ 

0.0049(13) 


82 


Pi P\P\ A A Pi ( <5 C\\ 

0.00449(39) 


84 


n o nT ! a -\ \ 

0.307(41) 


103 


n nincn/TO\ 

0.01050(73) 


Ne 


40 


n i/"*r"/in\ 

0.165(12) 


56 


n moi /oi7\ 

0.0181(37) 


63 


n nono / ^ /"> \ 

0.0393(46) 


83 


n pipipipio / 1 n\ 

0.00098(19) 


85 


n o n n / 7 \ 

0.800(65) 


104 


n nnoTo/7n\ 

0.00378(52) 


in n n n /"* o / n a \ 

19 0.0068(24) 


41 


n nn i n/To\ 

0.0910(73) 


57 


n nnTT/" i /^\ 

0.0077(16) 


64 


n i <r/i 1 \ 

0.145(11) 




Br 


86 


i nno/ on\ 

1.298(80) 




Tc 


n n n i n / i o \ 

20 0.126(13) 


42 


n n^i K / ^ o \ 

0.0415(43) 




Mn 


65 


n noo/i^?\ 

0.233(17) 


74 


n no^?T/o^?\ 

0.0377(87) 


87 


1.436(82) 


92 


n oni /nT\ 

0.301(37) 


Ol Pi OOC/00\ 

21 0.336(28) 


43 


n pit p t / o /i \ 

0.0167(24) 


50 


0.000280(81) 


66 


0.254(18) 


75 


Pi 01 T\ 

0.216(17) 


88 


1.197(65) 


93 


i o o / 1 n\ 

1.23(10) 


no n r"oo/^/i\ 

22 0.538(44) 




K 


51 


n nnfo/i -i \ 

0.0053(14) 


67 


n i to / i o \ 

0.173(13) 


76 


n ^ no / nn\ 

0.428(29) 


89 


n Tn/" 1 / /) o\ 

0.796(43) 


94 


n o t/ i 

2.87(16) 


23 0.247(45) 


39 


n nonn/o/" 1 \ 

0.0399(36) 


52 


n noo.i/no\ 

0.0384(93) 


68 


0.1044(85) 


77 


n i — i n / o t\ 

0.510(37) 


90 


n /i^n/n^\ 

0.440(24) 


95 


o Tn/nn\ 

3.79(20) 


Na 


40 


0.1328(92) 


53 


0.131(27) 


69 


0.0495(50) 


78 


0.387(34) 


91 


0.218(12) 


96 


3.51(19) 


22 0.0677(80) 


41 


0.198(14) 


54 


0.189(39) 


70 


0.0216(26) 


79 


0.347(25) 


92 


0.0978(53) 


97 


3.19(17) 


23 0.321(22) 


42 


0.162(12) 


55 


0.174(36) 


71 


0.00798(88) 


80 


0.223(14) 


93 


0.0446(25) 


98 


2.44(13) 


24 0.305(21) 


43 


0.1108(86) 


56 


0.098(20) 


72 


0.00311(38) 


81 


0.1365(83) 


94 


0.01625(97) 


99 


1.772(94) 


25 0.282(35) 


44 


0.0473(47) 


57 


0.058(12) 




Ga 


82 


0.0583(36) 


95 


0.00596(44) 


100 


1.115(61) 


Mg 


45 


0.0201(27) 


58 


0.0227(47) 


66 


0.096(12) 


83 


0.0244(16) 


96 


0.00228(28) 


101 


0.677(38) 


23 0.00559(46) 




Ca 


59 


0.0092(19) 


67 


0.238(18) 


84 


0.00833(60) 




Zr 


102 


0.351(21) 


24 0.1196(91) 


41 


0.0336(31) 


60 


0.00333(69) 


68 


0.288(20) 


85 


0.00221(23) 


86 


0.188(46) 


103 


0.165(11) 


25 0.271(27) 


42 


0.1245(89) 




Fe 


69 


0.245(17) 




Kr 


87 


1.060(81) 


104 


0.0705(58) 


26 0.394(44) 


43 


0.204(14) 


55 


0.116(10) 


70 


0.156(11) 


76 


0.0343(90) 


88 


1.51(10) 


105 


0.0287(33) 


27 0.242(72) 


44 


0.196(14) 


56 


0.226(17) 


71 


0.0830(69) 


77 


0.204(16) 


89 


1.87(11) 


106 


0.0101(23) 


Al 


45 


0.1149(87) 


57 


0.193(15) 


72 


0.0413(36) 


78 


0.484(33) 


90 


1.587(86) 






25 0.0015(38) 


46 


0.0560(51) 


58 


0.144(11) 


73 


0.0173(15) 


79 


0.574(44) 


91 


1.069(56) 
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TABLE IV: Isotopic cross sections for the production of elements ranging from Ru to Ba. The uncertainty includes both 
statistical and systematical errors. Where indicated in parenthesis, it applies to the last decimal digits. Underlined values are 
deduced from systematics. 



A a [mb] 


A 


a [mb] 


A 


a [mb] 


A 


n [ml>] 


.4 


a [mb] 


.4 


a [mb] 


Ru 


110 


0.698(38) 


118 


0.281(16) 


119 


5.85(30) 


118 


9.88(50) 


118 


0.121(18) 


95 1.174(85) 


111 


0.348(20) 


119 


0.1341(90) 


120 


4.15(21) 


119 


11.68(59) 


119 


0.395(23) 


96 2.99(17) 


112 


0.163(11) 


120 


0.0660(60) 


121 


2.79(14) 


120 


13.44(68) 


120 


0.954(52) 


97 4.42(23) 


113 


0.0744(59) 


121 


0.0288(46) 


122 


1.914(99) 


121 


14.39(73) 


121 


1.594(85) 


98 4.66(24) 


114 


0.0317(33) 


122 


0.0121(12) 


123 


1.232(64) 


122 


14.09(71) 


122 


2.86(15) 


99 4.29(22) 


115 


0.0140(84) 


123 


0.0050(10) 


124 


0.717(39) 


123 


13.88(70) 


123 


3.92(20) 


100 3.85(20) 




Ac 


124 


0.00194(37) 


125 


0.418(24) 


124 


12.44(63) 


124 


5.76(30) 


101 2.86(15) 


101 


0.480(53) 


125 


0.000690(61) 


126 


0.227(12) 


125 


11.00(56) 


125 


6.98(36) 


102 1.92(10) 


102 


1.95(12) 




In 


127 


0.106(11) 


126 


9.44(48) 


126 


9.10(46) 


103 1.142(62) 


103 


4.77(25) 


105 


0.092(24) 


128 


0.0517(59) 


127 


8.42(43) 


127 


10.27(52) 


104 0.632(36) 


104 


7.21(37) 


106 


0.760(59) 


129 


0.0198(11) 


128 


6.68(34) 


128 


12.71(64) 


105 0.321(19) 


105 


8.19(42) 


107 


2.72(15) 


130 


0.00560(33) 


129 


5.26(27) 


129 


13.63(69) 


106 0.157(11) 


106 


8.56(44) 


108 


5.81(31) 




Sb 


130 


4.26(22) 


130 


15.79(80) 


107 0.0581(53) 


107 


8.14(42) 


109 


8.49(44) 


110 


0.1330(87) 


131 


3.41(17) 


131 


17.34(88) 


108 0.0246(59) 


108 


6.80(35) 


110 


9.87(51) 


111 


0.664(38) 


132 


2.87(15) 


132 


21.2 ± 1.1 


109 0.0053(23) 


109 


5.40(28) 


111 


11.37(58) 


112 


2.23(12) 


133 


1.249(65) 


133 


23.3 ± 1.2 


Rh 


110 


3.72(19) 


112 


11.00(56) 


113 


4.85(25) 


134 


0.305(15) 


134 


31.5 ± 1.6 


97 0.961(85) 


111 


2.49(13) 


113 


10.14(52) 


114 


6.85(35) 




I 


135 


54.7 ±2.8 


98 2.69(16) 


112 


1.428(76) 


114 


8.43(43) 


115 


10.01(51) 


113 


0.00122(17) 




Cs 


99 4.91(26) 


113 


0.844(46) 


115 


6.46(33) 


116 


11.86(60) 


114 


0.0107(10) 


122 


0.145(11) 


100 5.63(29) 


114 


0.405(23) 


116 


4.79(25) 


117 


13.33(68) 


115 


0.0798(53) 


123 


0.275(18) 


101 5.75(30) 


115 


0.232(14) 


117 


3.24(17) 


118 


13.29(68) 


116 


0.448(26) 


124 


0.536(32) 


102 5.19(27) 


116 


0.1067(79) 


118 


2.05(11) 


119 


12.80(65) 


117 


1.247(67) 


125 


0.773(43) 


103 4.30(22) 


117 


0.0468(43) 


119 


1.331(70) 


120 


11.21(57) 


118 


2.31(12) 


126 


1.153(63) 


104 3.08(16) 


118 


0.0197(35) 


120 


0.755(40) 


121 


9.92(50) 


119 


4.09(21) 


127 


1.525(81) 


105 2.07(11) 


119 


0.0087(17) 


121 


0.438(24) 


122 


7.82(40) 


120 


6.00(31) 


128 


1.95(10) 


106 1.230(66) 


120 


0.00348(74) 


122 


0.249(15) 


123 


6.07(31) 


121 


8.48(43) 


129 


2.28(12) 


107 0.702(39) 


121 


0.001550(93) 


123 


0.1320(82) 


124 


4.76(24) 


122 


10.11(51) 


130 


2.81(15) 


108 0.336(20) 


122 


0.000420(37) 


124 


0.0580(61) 


125 


3.46(18) 


123 


12.76(65) 


131 


2.97(15) 


109 0.163(11) 




Cd 


125 


0.0290(41) 


126 


2.41(12) 


124 


13.24(67) 


132 


3.12(16) 




1U4 


1.26(11) 


izo 


n m oc/1 k\ 
U.U12o(10 ) 


1 97 

1 Z 1 




1 oc i 

1ZU 




LOO 


z. y*±i lo j 


111 0.0256(41) 


105 


3.80(22) 


127 


0.00420(23) 


128 


1.040(56) 


126 


14.97(76) 


134 


2.65(14) 


Pd 


106 


7.06(37) 




Sn 


129 


0.542(29) 


127 


16.51(84) 


135 


1.415(73) 


99 0.722(70) 


107 


8.49(44) 


108 


0.358(31) 


130 


0.290(15) 


128 


14.85(75) 


136 


0.499(28) 


100 2.49(15) 


108 


9.73(50) 


109 


1.628(94) 


131 


0.1414(75) 


129 


17.12(87) 




Ba 


101 4.94(26) 


109 


10.03(51) 


110 


4.34(23) 


132 


0.0396(20) 


130 


14.66(74) 


127 


0.0242(37) 


102 6.85(35) 


110 


9.04(46) 


111 


7.38(38) 




Te 


131 


16.38(83) 


128 


0.0326(47) 


103 6.91(36) 


111 


7.55(39) 


112 


9.48(49) 


111 


0.00777(63) 


132 


14.53(73) 


129 


0.0381(45) 


104 7.07(36) 


112 


5.79(30) 


113 


11.77(60) 


112 


0.0565(35) 


133 


15.40(78) 


130 


0.0267(39) 


105 6.11(31) 


113 


4.08(21) 


114 


12.49(64) 


113 


0.289(16) 


134 


21.9 ± 1.1 


131 


0.0228(36) 


106 4.79(25) 


114 


2.70(14) 


115 


12.43(63) 


114 


1.152(61) 


135 


23.6 ± 1.2 


132 


0.0147(29) 


107 3.35(17) 


115 


1.645(86) 


116 


11.45(58) 


115 


2.97(15) 




Xe 






108 2.22(12) 


116 


0.929(50) 


117 


9.58(49) 


116 


4.76(25) 


116 


0.00258(35) 






109 1.322(71) 


117 


0.547(30) 


118 


7.80(40) 


117 


7.37(38) 


117 


0.0188(19) 







two elements, cesium and barium, by charge exchange. 

Integrating all measured production cross sections, a 
value of 1393 ± 72 [mb] was obtained. Deducing the to- 
tal reaction cross section is not straight-forward, because 
the multiplicity of products from binary decays and mul- 
tifragmentation is larger than one. In addition, prod- 
ucts with Z < 3 are not observed. When we assume 
that products with Z < 3 are always accompanied by 
heaver fragments, the sum of the measured production 



cross section is an upper estimate of the total interaction 
cross section. The geometric cross section calculated by 
the model of Karol [4o| modified by Brohm [il|, which is 
1353 [mb], is consistent with our result. 

The production cross sections are summarised on the 
nuclide chart and represented with projections along the 
element and neutron number in Fig. [TTJ The general fea- 
ture is the U-shape of the element- and neutron-number 
distributions, which range over about two orders of mag- 
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FIG. 9: Production cross sections for the isotopes of elements ranging from Li to Ge. 



nitude. 



DISCUSSION 



A more detailed survey of the distribution of nuclide 
cross sections reveals the presence of staggering effects. 
In the region of the light nuclides, these are visible in the 
charge and neutron-number distributions (Fig. [II]) , and 
even more clearly in the distribution of production cross 
sections for specific values of TV — Z, as shown in Fig.[T3k: 
the staggering is also visible in the isotopic distribution 
of the heaviest residues like xenon and iodine, as shown 
in Fig. [H:. 



Without entering into the discussion on the phe- 
nomenology of the reaction, which is beyond the purpose 
of this report, we conclude by pointing out some physics 
cases which will be the subject of further research. 



A. Overview of the production of nuclides 

From the overview of the whole production of nu- 
clides we infer some general aspects of the decay pro- 
cess. The largest fraction of the reaction cross section re- 
sults in the production of heavy spallation residues, close 
to the mass of the projectile, which decays by evapora- 
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FIG. 10: Production cross sections for the isotopes of elements ranging from As to Ba. Cross sections indicated by circles were 
deduced from systematics. 



tion of mostly nucleons and clusters. This fraction pop- 
ulates the neutron-deficient side of the chart of the nu- 
clides, in agreement with the classic picture of spallation- 
evaporation: quite independently of the neutron enrich- 
ment of the hot spallation residues, the sequential evap- 
oration process tends to direct the decay towards the 
residue corridor in the region of heavy nuclides, 

this is located in the neutron-deficient side of the nuclide 
chart and corresponds to the situation which could be 
described as [47| dN/dZ = (Tn/Tz), where Tz and IV 
are the proton and neutron emission widths, respectively. 

Also, the evolution of the cross section with mass loss 
confirms the physical picture for the heavy spallation 



products drawn from previous experiments with neutron- 
rich systems [ll|. Fig.[l4]illustrates the close similarity of 
the 136 Xe( lj 4GeV) +P system with the heavy-residue pro- 
duction in 208 Pb(! a G e v) +p. The cross section evolves as 
a function of mass loss with very similar characteristics 
for the two systems: a steep fall describes the loss of the 
first few (less than ten) mass units; a plateau follows and 
changes gradually into an almost exponential descent be- 
yond around AA = 40 with an almost identical slope for 
both systems 208 Pb(iAGeV) +P an d 136 Xe(iAGeV) +p ■ 
Beyond A A = 70, abruptly, the two systems follow dif- 
ferent rules, according to the different decay processes 
which aliment the largest mass losses: symmetric binary 
splits are responsible for a significantly larger cross sec- 
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FIG. 11: (Color online) Nuclide production cross sections 
represented on the nuclide chart. Colours vary according to 
a logarithmic scale. The lower panel presents the projections 
of the residue-production cross sections along the atomic and 
neutron number. The charge-exchange contribution is also 
indicated in the projections. 



tion for the former system, which results in a large fis- 
sion hump, while, for the latter system, the intermediate- 
mass-fragment production models the mass distribution 
in the shape of a deep hollow. 

This sudden change of behaviour for large mass loss 
is also reflected in the average neutron-to-proton ratio 
(N/Z) of the residues, the evolution of which is shown in 
Fig. [321 In the region of the heaviest residues, above tin, 
the evolution with the element number of the quantity 
(N/Z), which indicates the ridge of the residue produc- 
tion, is compatible with an evaporation calculation where 
only neutrons, protons and alpha particles are emitted. 
As shown in the figure, the ridge calculated within this 
simplified evaporation pattern would evolve for light ele- 
ments towards the residue corridor. However, the quan- 
tity (N/Z) deviates from the neutron-deficient side of 



the chart of the nuclides around zirconium and migrates 
progressively towards the neutron-rich side for lighter 
residues. The lightest residues even populate the neutron 
rich side of the nuclide chart with respect to the stability 
valley. Values of (N/Z) connected with the neutron- rich 
side of the chart of the nuclides are a property of fission 
decays as well as the result of the multifragmentation of 
a neutron-rich source 14811. 



B. Staggering of nuclide cross sections 

The cross section oscillates between neighbouring nu- 
clides signing the presence of strong fine-structure ef- 
fects. As shown in Fig. [13k . the fine structure manifests 
over chains on nuclides with the same value of N — Z. 
An "even-odd" staggering characterises chains of nuclides 
having N — Z < 0. For N — Z > 0, the staggering is "even- 
odd", with a higher production of even nuclides, for all 
chains of even value of N — Z, but it reverses, with a 
higher production of odd nuclides, for odd values of N— Z 
for neutron-rich nuclides. The amplitude of the stagger- 
ing can be calculated with the procedure introduced by 
Tracy [49] , like in the analysis of the residue production 
from fragmentation reactions in Ref. [H(J. The result of 
the Tracy analysis applied to the cross sections of Fig. [13k 
is shown in Fig.fTgb. The highest amplitude of the "even- 
odd" staggering is measured for the chain N — Z = 0, and 
it may exceed 40%. A lower, but still large amplitude of 
the "even-odd" staggering characterises the other even 
value of N — Z: it evolves from 30% to 10% for the chain 
N — Z = 2, below chlorine. A similar evolution of the 
amplitude can be appreciated for the inverse "even-odd" 
staggering, which characterises the chains with odd val- 



1.5 
1.4 
1.3 
1.2 
1.1 
1, 



136 Xe| 








10 



20 



30 
Z 



40 



50 



60 
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ues of N — Z > 0, like N — Z = 3. In the regions of the 
chart of nuclides where the neutron-rich side is more pop- 
ulated than the neutron-deficient side, if fine-structure 
effects are present, they amplify in an "even-odd" stag- 
gering when the nuclide cross sections are projected on 
the neutron-number axis, and they compensate, resulting 
in a smoothed staggering, when projected on the atomic- 
number axis. The atomic- and neutron-number projec- 
tions shown in Fig. Qj] manifest this phenomenology: the 
strong "even-odd" staggering characterising the neutron- 
number distribution appears in correspondence with por- 
tions of the chart of nuclides where neutron-rich sides are 
most populated; these are two regions: the intermediate- 
mass fragments and the very heavy nuclides. These re- 
gions can be easily identified in Fig. [HI 

The staggering is understandable in relation with the 
decay process. Concerning the region of intermediate- 
mass fra gme nts, previous measurements of highly excited 
systems |50l . l5l| already revealed the same kind of stag- 
gering. It was suggested that the effect should be related 
to the reconstitution of fine-structure effects in the level 
densities during the cooling of the nucleus. In particu- 
lar, the formation of the intermediate-mass fragments is 
consistent with a process ending with at least few evap- 
oration steps which are responsible for the appearing of 
the staggering. Concerning the heaviest elements, mostly 
xenon and iodine, we observe an "even-odd" staggering 
along the isotopic distribution with an amplitude ap- 
proaching a value of 8% in the case of iodine, as shown 
in Fig. [T3h and Fig. [T3h . This feature is consistent with 
an almost complete survival of pairing correlations due 
to the low angular momentum which characterises the 
production of the heaviest fragments [H3.l53j|. 



C. Hints on the reaction mechanism 

The kinetic energies of the fragments give hints on the 
reaction kinematics. As shown in Fig. the part of the 
experimental spectrum related to the intermediate-mass- 
fragment production deviates in favour of very high ki- 
netic energies from the systematics, which describes the 
contribution from the heaviest residue of one decay se- 
quence, characterised by a larger number of small recoils 
from evaporation of particles and light fragments. In the 
part of the spectrum where the mass of the residues is 
lower than half the mass of the projectile, the deviation 
from the systematics is due to the predominant contribu- 
tion of the light partner of a binary decay or a product 
of a multifagmentation event characterised by a high ki- 
netic energy. We may remark that also the U-shape of the 
mass distribution, like the extension of the production to 
rather neutron-rich nuclides, are properties of multifrag- 
mentation as well as fission below the Businaro-Gallone 
point. 

Concerning the evolution of production cross section 
as a function of the mass loss, shown in Fig. HH we pro- 
pose a comparison of the slope of the mass-loss distribu- 
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FIG. 13: (Color online) (a). An "even-odd" staggering in the 
production cross section of the residues manifests for chains of 
even value of N — Z with a higher production of even nuclides, 
and it reverses, with a higher production of odd nuclides, for 
odd values oi N — Z for neutron-rich nuclides, (b). Stagger- 
ing amplitude of the cross sections shown in (a) analysed by 
the method of Tracy [49l |. It is positive when the staggering 
is "even-odd" and negative when the staggering reverses, (c). 
"Even-odd" staggering in the production cross sections of the 
isotopes of xenon and iodine in 136 Xe(i a Gev) +P • (d). Stag- 
gering amplitude of the cross sections shown in (c) analysed 
by the method of Tracy [4j|. 



tion for the heavy-fragment side of the minimum, with 
a simplified evaporation calculation limited to the emis- 
sion of neutrons, protons and alpha particles. This side 
of the spectrum may be taken as a reference for test- 
ing the validity of cascade models, due to its sensitivity 
to the excitation energy introduced into the system by 
the collision. However, in this energy range, the function 
we estimated (dashed line in Fig. [T4|) underpredicts the 
measured spectrum. The high kinetic energies of the light 
residues and the U-shape of the spectrum, suggest that 
the intermediate-mass-fragment production is mostly al- 
imented by asymmetric splits of heavy nuclei and, there- 
fore, the difference between the measured descent of the 
mass-loss distribution in the range AA = 40 to 70 and 
the simplified evaporation calculation corresponds to the 
heavy partners of light fragments originating from asym- 
metric splits. 

It is possible that this explanation, already pro- 
posed for lighter systems like 56 Fe(iAGeV) + P [3 
is rather general and describes also heavier systems 
at the energy range of 1 GeV per nucleon like, for 
instance, the 208 Pb(iAGeV) + P system. The speci- 
ficity of 136 Xe( lj 4G e v) + p as compared to the previ- 
ous study of lead consists in the possibility to anal- 
yse the intermediate-mass-fragment production. The 
comparison of the two systems suggests that a simi- 
lar intermediate-mass-fragment production could charac- 



1G 



terise also the Pb(i ^GeV) +P system and increase the 
production cross section in the range A A = 40 to 70. 

VI. CONCLUSIONS 

The production cross sections of fully identified nu- 
clides produced in the system 136 Xe(i a GeV) +P were mea- 
sured covering the element range from lithium to barium. 
The measurement of the kinetic energies imparted to the 
emitted residues completes the set of experimental data 
which is of high relevance both for applications, like the 
damage and activation of irradiated materials, and, more 
generally, for the modelling of spallation reactions. 

The present measurement enables the first full sur- 
vey of intermediate-mass-fragment cross sections mea- 
sured together with the heavy-fragment production in 
the decay of a neutron-rich system in the fissility region 
where asymmetric splits are predominant. The process of 
heavy-fragment formation in the system 136 Xe(i 7 iGeV) + 
p, determined mostly by sequential evaporation presents 
a phenomenology already observed in the decay of 
neutron-rich systems and establishes a direct connection 
with the system 208 Pb(i J 4GeV) + P > measured in a former 
experiment [aV The intermediate-mass-fragment produc- 
tion recalls the phenomenology observed in the system 



56 Fe(i J 4GeV)+P > resumed by the U-shape of the mass dis- 
tribution and the very high kinetic energies which char- 
acterise the process [la]. These comparisons even sug- 
gest extending the description of the intermediate-mass- 
fragment formation, firstly proposed for light system and 
confirmed in the present work for the decay of xenon, also 
to the heavier neutron-rich systems, in the fissility region 
where symmetric splits are largely more probable than 
asymmetric splits. As an extreme case, the light nuclide 
formation characterising 238 U(i J 4GeV) + P UM ma y be 
explained by a similar scenario. 



In this respect, the system Xe^^QeV) +P is a rep- 
resentative system for studying the intermediate-mass- 
fragment formation in connection with the spallation pro- 
cess induced by protons in the 1 A GeV energy range on 
neutron- rich targets; we propose it as a benchmark for 
testing the validity of models for spallation reactions over 
a large range of residue masses and charges. A correct 
identification of the process of formation of intermediate- 
mass fragments requires a dedicated analysis of the re- 
action kinematics and of the role of the Coulomb repul- 
sion, as well as an exclusive measurement of correlation 
observables, and is a perspective for future research. 
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